Exploring the phase space of the quantum delta kicked accelerator 
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We experimentally explore the underlying pseudo-classical phase space structure of the quantum 
delta kicked accelerator. This was achieved by exposing a Bose-Einstein condensate to the spatially 
corrugated potential created by pulses of an off-resonant standing light wave. For the first time 
quantum accelerator modes were realized in such a system. By utilizing the narrow momentum 
distribution of the condensate we were able to observe the discrete momentum state structure of a 
quantum accelerator mode and also to directly measure the size of the structures in the phase space. 
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For more than a century the study of chaotic phenom- 
ena has been recognized as being crucial to developing a 
fuller understanding of nature. One aspect of this study 
which was missing until relatively recently was experi- 
mental scrutiny of quantum systems which in the classical 
limit exhibit chaotic behavior. Theoretical work on this 
subject had largely concentrated on the investigation of 
idealized systems such as the quantum delta-kicked rotor 
(QDKR) which were already well known from extensive 
work in the classical regime jlj. The experimental study 
of this system gained new impetus through its realiza- 
tion using laser cooled atoms exposed to a corrugated 
otential from a pulsed off-resonant standing light wave 
This system has subsequently led to many discover- 
ies in the field of quantum chaos including observation of 
quantum resonances [3,J3|, dynamical localization 0, 0| 
and quantum diffusion [a, Lfll ■ In the theoretical descrip- 
tion of the QDKR the effective value of Planck's constant 
scales with the time between the pulses Therefore, 
to achieve classical correspondence in which S — > 0, the 
time between pulses needs to be close to zero. However, 
for technical reasons this value can not be made arbitrar- 
ily small in experiments. 
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FIG. 1: Experimental momentum distributions showing a 
QAM in a BEG exposed to a series of kicks from a stand- 
ing light wave in the free falling frame. The kicking period 
was 72.2 fis and the momentum distributions are displaced as 
a function of the number of kicks. The accelerator mode is the 
collection of momentum states that appear to move towards 
the upper right. The inset shows the momentum distribution 
for a pulse number of n=21. 



Recently, Fishman, Guarneri, and Rebuzzini (FGR) 
|ir)l | have shown that this difficulty can be circumvented 
for kicking periods close to a quantum resonance time. 
These resonance times occur at integer multi- 

ples of the Talbot time, a time interval during which 
plane waves with certain equally spaced momenta in the 
kicking direction can acquire a phase factor which is an 
integer multiple of 2tt. This is analogous to the Talbot 
effect in optics j^^. Using this method, FGR have stud- 
ied the quantum delta-kicked accelerator (QDKA) which 
can be created by taking the QDKR and adding a lin- 
ear potential along the direction of the standing wave. 
They showed that the effective value of Planck's con- 
stant scales with the deviation of the pulse separation 
time from certain non-zero resonance times, thus mak- 
ing it feasible to make the effective Planck's constant 
very small. Therefore, if the time between pulses is cho- 
sen close to these resonance times, a pseudo-classical ap- 
proach can be adopted to study the system. 



Perhaps the simplest way of experimentally realizing 
the QDKA is by applying the pulsed standing wave in 
the direction of gravity 0, . This experiment has led 
to the discovery of quantum accelerator modes (QAMs) 
near the resonance times. One of the most important 
characteristics of QAMs is that they are comprised of 
atoms which show a linear momentum growth with pulse 
number in a freely falling frame ffsj . It has been shown 
that QAMs are quantum nondissipative counterparts of 
mode locking [l6j. They have also been suggested for 
use in the preparation of well defined initial conditions 
for quantum chaos experiments FGR attributed 

the QAMs to the existence of stability islands in the 
pseudo-classical phase space. These studies have shown 
that this underlying phase space has a complex structure 
which is highly sensitive to the experimental parameters. 
However, the broad momentum distribution of the laser 
cooled atoms which have been used so far to study this 
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system have prevented the examination of the local struc- 
tures in the phase space. 

In this letter we report on the realization of QAMs 
using a Bose-Einstein condensate (BEC) of rubidium 87 
atoms and the exploration of the pseudo-classical phase 
space structure of the QDKA. Figure^shows experimen- 
tally observed momentum distributions as a function of 
the number of standing wave pulses applied to a BEC. 
This figure demonstrates that the QAM gains momentum 
linearly as the number of pulses increases. Note that this 
is the first time that it has been possible to determine 
that the QAM is made up of several distinct momentum 
states as originally postulated in Ref. This quanti- 
zation of momentum is observable because the initial mo- 
mentum uncertainty of the condensate was much smaller 
than two photon recoils. 

The Hamiltonian of the delta kicked accelerator for an 
atom with mass m can be written as, H — /2m + 
mg'x + -\- cos(Ga;)] ^„ 5{t — nT), where p is the 

atomic momentum along the standing wave, g' is the 
component of the gravitational acceleration in the direc- 
tion of the standing wave, G = 2k the grating vector 
where k is the light wave vector, n the pulse number, T 
the pulse period, and Umax is the well depth of the stand- 
ing wave due to the light shift. The net effect of the time 
dependent potential is to distribute the condensate into 
different momentum states via a diffractive process. The 
population of these momentum states is determined by 
\Jn'i<Pd)\'^, where Jn'{4>d) is an n'th order Bessel function 
of the first kind with argument (/>d — Umax^t / {2h) , the 
phase modulation depth [lj| . It is useful to consider this 
system for pulse periods close to integer multiples of the 
half Talbot time, Te ^ i x 2Tim/hG'^(^ I x 33.3^s for 
Rb87 atoms). With this restriction the system can be 
described by the classical map [loj . 
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where e = 271 £(T/Ti — 1) is a small number, £ is any pos- 
itive integer number, and K = \e\(l)d- The dimensionless 
J and 9 parameters are defined as, 

e = G X mod(27r) 
Jn = In+sgn{e)[7Ti + T{P + m^ + T]/2)] (2) 

where p/{hG) = + /3, /3 is the fractional momentum, 
T = hTG'^/m, and r/ = mg'T/{hG). Figure shows a 
typical phase space portrait for the map of Eq. with 
(j)d = 1.4, T = 29.5^s and e — —0.72. Perhaps the most 
important feature of this plot is the existence of a stable 
fixed point surrounded by an island of stability. If the size 
of these islands is large enough to capture a significant 
fraction of the wavepacket they give rise to observable 
accelerator modes. According to this model the momen- 
tum gain of an atom in a period p accelerator mode after 
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FIG. 2: Phase space unit cell for the QDKA map of Eq. Q 
with (jjd = 1.4 and T = 29.5 /xs. A stable fixed point with 
(p,m) = (1,0) exists at J = and 61 = 0.0887. This stable 
fixed point is surrounded by the stability island in which the 
quantum accelerator mode will be created if the atomic initial 
conditions are inside the island. The momentum width of 
the condensate, AJbec and the stability island AJisiand, are 
shown with the solid white and black lines respectively. 



n kicks is given by. 
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where m is an integer and (p,m) specifies a particular 
accelerator mode [llj, llal . 

The initial conditions of atoms that can be accelerated 
is another important property of the accelerator modes. 
The limited range of these conditions is a consequence of 
the fact that the stability islands do not cover the whole 
unit cell of the phase space. Furthermore, the initial 
momentum required for an accelerator mode to emerge 
is periodic. This can be seen by using Eq. (O and the 
fact that J has a periodicity of 27r. This is equivalent 
to a momentum periodicity of Ap — 2-KhG/T. Observ- 
ing this phase space structure requires that the atomic 
momentum distribution be narrower than 2ttHG/t. This 
implies a temperature of w 450 nK for Rb87 atoms ex- 
posed to a pulsed standing wave of 390 nm wavelength 
and pulse period close to the Talbot time {i = 2). To this 
extent, a BEC is an ideal candidate (lij. The momen- 
tum width of a condensate is shown with two white lines 
in Fig. 12 Note that this width is smaller than the mo- 
mentum extent of the island, indicating that the momen- 
tum resolution of the experiment is more than sufficient 
to clearly detect and identify a stability island from the 
chaotic background. However, for experiments utilizing 
cold thermal atomic samples, the momentum distribu- 
tion is significantly wider than HG. Although this wide 
momentum distribution makes it relatively easy to ob- 
serve the accelerator modes, there is no direct way of 



3 



examining the structure of the phase space. 

In order to use BEC to explore the phase space of 
the kicked accelerator we used the experimental setup 
described in detail by Ref . 0| . A standard six-beam 
MOT was used to trap about 50 x 10^ atoms which were 
loaded into the trapping potential created by a focused 
CO2 laser beam. In order to optimize the loading effi- 
ciency, the waist of the beam was chosen to be 100 ^m 
[2H . Typically 4 x 10^ atoms were trapped at the focus 
of the CO2 laser beam. Subsequently, one of the lenses 
through which the beam passed was moved 17 mm in 
1 s so as to reduce the beam waist to 12 /zm and com- 
press the trap for efficient evaporative cooling j^^- The 
power in the beam was then reduced to 100 mW in 4 s 
to create a pure condensate with 50000 atoms in the 
5S1/2F = l,mp = state. The CO2 laser was then 
turned off and after a variable time interval the kicking 
potential was turned on. Varying this time allowed the 
BEC to fall under the influence of gravity, thus changing 
the momentum of the condensate at the commencement 
of the kicks. The light for the kicks was obtained by 
passing the light from the MOT laser through a 40 MHz 
acousto-optic modulator (AOM). This light was 6.7 GHz 
red detuned from the 5S1/2F = 1 5P3/2F = 2 transi- 
tion of Rb87 and propagated at 41° relative to the verti- 
cal direction. The beam size was 1 mm, such that (pd did 
not change appreciably while the BEC was interacting 
with the series of kicks. The phase modulation depth, 
(f>d, was inferred by comparing the population of the first 
and zeroth order momentum states after one pulse. The 
temporal profile of the standing light wave was controlled 
by periodically switching the AOM on in order to create a 
sequence of pulses each 250ns in length. The momentum 
distribution was measured using a time of flight method. 
That is, the condensate expanded for a controlled time 
interval, typically 9 ms, and was then destructively im- 
aged using an absorptive technique. 

To observe the pseudo-classical phase space structure 
of the QAMs, a series of data were taken for pulse periods 
near both Talbot and half- Talbot times. Figure^ shows 
a typical data set taken at (a) T = 61^s, (b) T = 72.2/xs, 
(c) T = 28.5^s, and (d) T — 37.1^s pulse separations for 
different values of the BEC's initial momentum, (a) and 
(b) occur on either side of the Talbot time at T = 66.6/iS 
while (c) and (d) occur on either side of the half- Talbot 
time at T = 33.3/iS. The data in Figure|2|was created by 
horizontally stacking 60 time-of-flight images of the con- 
densate, each for a different initial momentum. These 
data conffim the periodicity of the QAMs with momen- 
tum. Furthermore, the data of Fig.O provides a direct 
way to validate the theoretical prediction of the island 
size. To do so, the data of Fig.^was summed along the 
initial momentum axis. Apisiand was then determined by 
measuring when the accumulated signal of the QAM had 
dropped to 1/e of its maximum value. The theoretical 
values were inferred by plotting the map of Eq. Q for 
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FIG. 3: Experimental momentum distributions showing the 
sensitivity of the QAMs to the initial conditions. The data 
is taken after applying 10 kicks with pulse intervals of (a) 61 
fis and (b) 72.2 and 30 kicks with pulse intervals of (c) 
28.5 /is and (d) 37.1 /is. The larger number of kicks at the 
half- Talbot time was necessary due to the smaller momentum 
transfer per pulse. The initial momentum was changed by 
applying the kicking potential with variable time delays after 
releasing the condensate from the CO2 laser. The slope of the 
data seen in this figure is caused by the momentum gain due 
to gravity. The doted lines denote the position of the QAM. 
Note that Apisiand is related to the size of the stability islands 

by Apisland = fiGAJisland/T- 



the corresponding experimental values of K. The exper- 
imental and theoretical values for the momentum extent 
of the islands are given in Fig.^ near (a) half- Talbot and 
(b) Talbot times. The circle and asterisk signs are the ex- 
perimental and theoretical values for A Jisiand (as defined 
in Figs. m and O. It can be seen that the experimental 
values are very close to the theoretical predictions. Note 
that for values of iiT > 2 at half- Talbot time, the stability 
island elongates in J and becomes narrow in 6. This be- 
havior reduces the effective overlap between the BEC's 
wavefunction and the stability island and consequently 
the QAMs were not visible in Fig.Qfb) for higher values 
of K. Figure|2| also shows that there can be little overlap 
between the initial conditions that will populate a QAM 
at two different values of T. This behavior particularly 
affects what happens in experiments in which the mo- 
mentum distribution is measured as the pulse period is 
scanned across a resonance time. Unlike the experiments 
with cold atomic samples where the QAMs on both sides 
of a resonance could be populated [l^ , in the case of the 
condensate, only the QAMs which have significant over- 
lap with the condensate wavefunction will be observable. 
This can be seen in Fig.El where we performed a scan 
of pulse period across the Talbot time for two different 
initial momenta. The initial momentum for Fig.[SJa) was 
set to 1.2hG such that the QAMs were efficiently loaded 
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FIG. 4: Experimental data showing the momentum range in 
which the QAMs appear and the corresponding theoretical 
predictions, (a) shows the experimental (circle) and theoret- 
ical (asterisk) values for AJiaiand near half- Talbot time and 
(b) show the same quantities near Talbot time. Figure El il- 
lustrates the method used to experimentally infer AJisiand- 
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FIG. 5: Experimental momentum distributions showing con- 
trollable loading of the (p,rn) =(1,0) QAM. The data was 
taken by applying 10 kicks for a range of kicking pulse peri- 
ods. In (a) the initial momentum was set to 1.2hG so as to 
efficiently create the QAM at pulse periods greater than the 
Talbot time while in (b) it was set to l.bhG to efficiently cre- 
ate the QAM at pulse periods smaller than the Talbot time. 
The curves show the position of the QAM calculated with 
Eq. (3). 



at pulse periods near T — 72/is, whereas in Fig.|SJb), the 
initial momentum was set at l.SftG to mainly populate 
the QAMs around pulse periods near T — 61/is. As can 
be seen, the QAM with indices (p,tn) =(1,0) appear at 
pulse periods greater than the Talbot time in Fig.|3fa), 
whereas in Fig.jSJb) the (1,0) QAM mostly appears at 
pulse periods smaller than the Talbot time. Note that 
this is the first time that it has been possible to selec- 
tively populate an island at a particular position in phase 



space. 

In conclusion our experiments have demonstrated the 
feasibility of observing quantum accelerator modes using 
a BEC. Using a BEC we were able to examine the under- 
lying pseudo-classical phase space structure of the quan- 
tum delta kicked accelerator. These experiments pave the 
way towards further investigation of more complex sys- 
tems like the kicked harmonic oscillator [2^ and dy- 
namical tunneling ^,, ,26] that have received little study. 
This experiment also opens the door to experimental ob- 
servation of many phenomena related to quantum chaos. 
For example, the high momentum resolution of the exper- 
iment could bring the possibility of observing bifurcation 
of the stability islands to a practical level. This should 
lead to a better understanding of the transition to chaos 
in a quantum system. 



[1] 

[2] 
[3] 
[4] 
[5] 

[6: 

[7 

[s: 

[9 
[10 

[11 
[12 
[13 

[14 
[15 
[16 
[17 

[18 

[19 
[20 

[21 

[22; 

[23; 

[24; 

[25; 

[26; 



A.J. Lichtenberg and M.A. Lieberman, "Regular and 
chaotic Dynamics" Springer- Verlag, New York (1992). 
F.L.Moore etal., Phys. Rev. Lett. 75, 4598 (1995). 
C.Ryu etal., Phys. Rev. Lett. 96, 160403 (2006). 

F. L.Moore etal., Phys. Rev. Lett. 73, 2974 (1994). 

J. Ringot, P. Szriftgiser, J.C. Garreau, and D.Delande 
Phys. Rev. Lett. 85, 2741 (2000). 

H. Ammann, R. Gray, I. Shvarchuck, and N. Christensen, 
Phys. Rev. Lett. 80, 4111 (1998). 

G. J.Duffy etal., Phys. Rev. E 70, 056206 (2004). 
W.H.Oskay, D.A.Steck, V.Milner, E.G. Klappuf, 
M.G. Raizen, Opt. Comm. 179, 137 (2000). 

M. Sadgrove, S. Wimberger, S. Parkins, and 
R. Leonhardt, Phys. Rev. Lett. 94, 174103 (2005) 
S. Fishman, I. Guarneri, and L. Rebuzzini, J. Stat. Phys. 
110, 911 (2003); ibid, Phys. Rev. Lett. 89, 084101 (2002). 

C. F. Bharucha etal., Phys. Rev. E 60, 3881 (1999). 
M.B.d'Arcy etal., Phys. Rev. Lett. 87, 074102 (2001). 
M.V. Berry and E. Bodenschatz J. Mod. Opt. 46, 349 
(1999). 

R.M. Godun etal., Phys. Rev. A 62, 013411 (2000). 

M.K. Oberthaler etal., Phys. Rev. Lett. 83, 4447 (1999). 

A. Buchleitner etal., Phys. Rev. Lett. 96, 164101 (2006). 

M.B.d'Arcy, R.M. Godun, D. Cassettari, and 

G.S.Summy, Phys. Rev. A 67, 023605 (2003). 

S. Schlunk, M.B. d'Arcy, S.A. Gardiner, and G.S. Summy, 

Phys. Rev. Lett. 90, 124102 (2003). 

L.Deng etal., Phys. Rev. Lett. 83, 5407 (1999). 

P. Ahmadi, G. Behin-Aein, B.P. Timmons, and 

G.S.Summy, J. Phys. B 39, 1159 (2006). 

P. Ahmadi, B.P. Timmons, and G.S. Summy, 

Phys. Rev. A 72, 023411 (2005). 

T. Kinoshita, T. Wenger, and D.S.Weiss, Phys. Rev. A 
71, 011602(R) (2005). 

S.A. Gardiner, D. Jaksch, R. Dum, J.I. Cirac, and 
P. ZoUer, Phys. Rev. A 62, 023612 (2000). 
G.J.Duffy, A.S.Mellish, K.J. Challis, and A.C.Wilson, 
Phys. Rev. A 70, 041602(R) (2004). 

D. A.Steck, W.H.Oskay, and M.G. Raizen, Science 293, 
274 (2001). 

W.K. Hensinger et al. Nature 412, 52 (2001). 



